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X 1 Abstract 

We calculate the semileptonic branching ratio and the charm counting of inclusive B 
meson decays using the perturbative QCD formalism. For the nonleptonic decays, we 
employ the modified factorization theorem, in which Wilson coefficients evolve with the 
characteristic scales of the decay modes. It is found that the decay rate of the single-charm 
mode b — > cud is enhanced, and a lower B$l is obtained without increasing n c . We predict 
a larger b — > ctv branching ratio compared to that from the conventional heavy-quark- 
effective-theory based operator product expansion. Our result of the B meson lifetime is 
also consistent with the data. 
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With rapid theoretical and experimental developments, the heavy flavor physics has 
reached the level of accuracy that the data can be used to test the standard model. Es- 
pecially, the various quantities related to the heavy meson decays have been calculated 
and compared with experimental data. In spite of many successes, a discrepancy exists 
consistently between the measurements and the theoretical predictions of the inclusive 
semileptonic branching ratio I?sl and the charm counting n c of the B meson decays: The 
latest results from CLEO, measured at the T(4s) resonance, are B SL = (10.19 ± 0.37)% 
and n c = 1.12 ± 0.05 |], while LEP reports, at the Z resonance, B SL = (11.12 ± 0.20)% 
and n c = 1.20 ± 0.07 |2|]. The naive parton model, which coincides with the leading-order 
heavy-quark-effective-field-theory (HQEFT) calculation, gives -Bsl ~ 13% HJ. Convention- 
ally, the inclusive B meson decays are described by a systematic HQEFT based Euclidean 
space operator product expansion (OPE) of relevant hadronic matrix elements in the in- 
verse power of the b quark mass mj. The assumption of global and local quark-hadron 
dualities is employed to justify this approach. The l/m| corrections (there are no l/ra& 
corrections) have been investigated and found to be less than 5% of the leading-order pre- 
dictions Q. Hence, it was argued that the values of £>sl < 12.5% can not be accommodated 
by QCD theory ||, and this discrepancy could be an indication of "new physics" beyond 
the standard model. For example, an enhanced b —>■ sg vertex has been proposed to explain 
the anomaly |J. 

It is quite possible that the assumption of the local duality fails at the l/m& level. Re- 
cently, there was theoretical indication that l/m& corrections might appear in nonleptonic 
decays [0. Altarelli et al. has observed that if m& in the phase space factor was replaced 
by the B meson mass Mb, and by M\ b in the decays, due to the breakdown of the 
local duality, the predictions agree well with the data of -Bsl and of the lifetime ratio ||. 
This naive modification increases the nonleptonic decay rates, and thus reduces -Bsl- On 
the other hand, though the nonperturbative corrections are small, it was pointed out that 
perturbative corrections may be significant ||. Bagan et al. have computed the 0(a s ) 
corrections with the effects of the charm mass included ||, and found that -Bsl is indeed 
lowered with the increase of the b — > ccs rate. However, the predictions are sensitive to 
the choice of the renormalization scale fi [111: .Bsl = 12.0 for fi = m&, and I?sl — 10.9 
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for fi = rrib/2. Furthermore, an enhanced b — > ccs mode, while lowering -Bsl, enlarges n c 
to around 1.20, which, though consistent with the LEP data, becomes significantly higher 
than the CLEO data. This is sometimes called the "missing charm puzzle" . 

In this paper we shall analyze the inclusive B meson decays based on perturbative QCD 
(PQCD) factorization theorems fnf . It will be shown that the experimental data of Bsl 
and n c can be accommodated simultaneously within the framework of the standard model, 
without resort to new physics. In a series of recent works jl2, [L| we have developed a PQCD 
formalism for the inclusive semileptonic decay B — > X u ev, whose decay rate is written as the 
convolution of a hard subamplitude with the B meson distribution function and the w-quark 
jet function. The idea is briefly explained as follows: Loop corrections produce infrared 
divergences, which signal the sensitivity to long-distance physics. The soft divergences 
associated with vanishing loop momenta survive after summing over all diagrams of the 
same order of a s . These divergences are absorbed into the B meson distribution function, 
which describes the probability that the b quark carries some fraction of the B meson 
momentum. This function has been identified as the soft function associated with the 
Fermi motion of the heavy quark inside the heavy hadron, and as the outcome of the 
resummation of nonperturbative power corrections in HQEFT [|14j] . Near the end points of 
the decay spectra, the invariant mass of the u quark vanishes, and the collinear divergences 
due to radiative gluons parallel to the u quark become important. Large double logarithms 
from the overlap of the soft and collinear divergences then arise, and are absorbed into 
the jet function. The remaining finite piece of the radiative corrections is grouped into the 
hard subamplitude. 

The hard subamplitude, dominated by short- distance contributions, is evaluated using 
perturbation theory. The end-point double logarithms contained in the jet function are sys- 
tematically resumed into a Sudakov factor. It turns out that the Sudakov factor suppresses 
the long-distance contributions, and renders the perturbative calculation more reliable and 
self-consistent. The nonperturbative B meson distribution function is universal, and de- 
termined by other experiment data, such as the spectrum of the inclusive radiative decay 
B — > X s 7 |15| . Note that the analysis of inclusive heavy meson decays suffers the infrared 
renormalon ambiguity, which appears as power corrections starting with 0(1 /ml). This 
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ambiguity can be parametrized as an exponential factor associated with the distribution 



function, which has also been fixed from the best fit to the B — > X s "y spectrum [15]. The 
above prescriptions specify clearly the treatment of each subprocess in the formalism. The 
advantage of our approach lies in that it is formulated according to the meson kinematics, 
instead of the parton kinematics as used in HQEFT. Hence, the replacement of mj by Mb 
postulated in M is implemented naturally. The kinematic gap between the parton picture 
and the hadron picture is also removed. 

We apply the above PQCD factorization theorem to the charmed decays, 

B(P B ) ^X c + l( Pl ) + v(p„), (1) 

which dominate the inclusive semileptonic B meson decays. The momenta are expressed, 
in terms of light-cone coordinates, as 

(M B M b 



B 



0± , p/ = (Pi,Pi ,0_l), p v = (pt,P v ,Pv±), (2) 



where the minus component p~[ vanishes for massless leptons. We choose Ei, q 2 and qo, or 
equivalently, the three dimensionless quantities, 



_ 2E l _ q 2 _ 2q 

X ~W B ' y -M%' y ° ~ ~Mb ' 



as the independent variables, E\ being the lepton energy and q = p\ + p v the momentum 
of the lepton pair. Their kinematic ranges are 

2m; 

_u x <i_^ + « ; 

M B 

x + \Jx 2 — Act 

en < y < a + (1 + a — p — x) 



2 — x — \/x 2 — 4a ' 



2(y — a) 

x + i ^ yo < 1 + y-fr ( 4 ) 

x + V£ — 4a 



with the constants 



™i a M D 
a = — V, P = —jtr , (5) 

M 2 B M 2 B 



and the lepton mass m; and the D meson mass Mp. Mp is introduced as the minimal 
invariant mass of the decay product X c . It is easy to check that the above ranges reduce 



to those for the decay B — ► X u ev as mi, Mb — > fL2 



In the parton picture the B meson is composed of a b quark with momentum p& = 
Pb — p, where p = (p + , 0, pj_) is the momentum carried by other light degrees of freedom. 
Since the c quark is massive, there are no collinear and thus no large double logarithms 
associated with it, differing from the decay B — > X u ev fll2 |. That is, QCD corrections 
to the charmed semileptonic decays are weaker. A convenient factorization scale is the 
transverse momentum p±, or its conjugate variable b, of the b quark inside the B meson 
T6fl . 1/b marks the scale separation, below which physics is absorbed into the B meson 



distribution function /, and above which physics is absorbed into the hard b quark decay 
subamplitude H. Therefore, we write the decay rates as, 

1 d 3 T M 2 B n 



H 0) dxdydy 



dz \ dbbf(z,b)J c (b)H(z) 



(6) 



with = (G F /lQTT 3 )\V c b\ 2 M B . The momentum fraction z is defined as z = P^/Pb 
I — p + /Pg. z is equal to the maximum 1 as the b quark takes all the B meson momentum 
in the plus direction. The minimum of z is determined by the condition p 2 > m 2 , p c and 
m c being the c quark momentum and mass, respectively, which leads to 
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(7) 



The function J c denotes the final-state cut, which is in fact part of the hard subamplitude. 
We factorize it out in order to compare Eq. (H) with the formulas in the nonleptonic case 
below. J c and H in momentum space are given by 



Jr. 



Pc 



5 M 2 
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The universal distribution function determined from the decay spectrum of B — > X s 7 is 



written as 15 



0.02647^(1 - z 



>2 



6 ) = [(, 0.95)* ; 0.0034^ ex P (-0.0125M^-) , (10) 
where the exponential parametrizes the infrared renormalon contribution at large 6. 

We then extend our formalsim to the nonleptonic decays B — > X c5 and B — > X c , 
which are dominated by the 6 quark decays b — > ccs and 6 — > cud, respectively. For 
simplicity, we ignore the pretty small charmless processes. To simplify the calculation, 
we route the transverse momentum p± of the b quark through the c quark jet as in the 
semileptonic cases, such that the c quark momentum is the same as before. For kinematics, 
we make the correspondence with the c (u) quark carrying the momentum of the massive 
(light) lepton r (e and /i) and with the s and d quarks carrying the momentum of v. 
The dimensionless variables are then defined exactly as Eq. (]3|). The nonleptonic decays 
involve more complicated strong interactions characterized by three scales: the W boson 
mass Mw, the characteristic scale t of the B meson decay, and the factorization scale 
1/6. The hard gluon exchanges among the quarks generate logarithms of Mw, and the 
associated physics is factorized into a harder function [pjj . Collinear divergences and soft 
divergences from radiative corrections exist simultaneously, when the light quarks, such as 
u, d, and s, become energetic. The resultant double logarithms ln 2 (p6), p being the jet 
momentum defined later, are factorized into the corresponding jet functions. This is how 
the additional scales M\y and p are introduced into the nonleptonic decays. There are also 
soft gluon exchanges among the quarks, whose effects are negligible in bottom decays [ |18|] . 
Similarly, physics with the scale below 1/6 is absorbed into the (soft) distribution function. 

We employ the three-scale factorization theorem developed in [|I7| , that embodies both 
the PQCD factorization theorem stated above and effective field theory. According to the 
above explanation, the decay rates are expressed as the convolution of the harder, hard, 
jet and soft subprocesses, 

r r r d 2 b ~ 

— = Ml J dxdydy dz J -^—^H r {M w , p)H(z, t, p) J[ J t (6) JJ Jj(b,Pj, n)S(z, b,fi) , (11) 

where the renormalization scale /i denotes the inclusion of QCD corrections. The functions 
Ji, i — c and c for the 6 — > ccs mode and i = c for the 6 — > cud mode, do not contain double 
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logarithms, and represent the final-state cuts. The light-quark jet functions Jj, j = s for 
the b — > ccs mode and j = u, d for the b — *> cud mode, contain double logarithm. These 



double logarithms are resummed into a Sudakov factor [12, 13 



Jj(b,pj,n) = eocp[-2s(p i ,6)]J i (6,/i) 



(12) 



where the upper bound of the Sudakov evolution is chosen as the sum of the longitudinal 
components of pj, 

Pj = (Pj + Pj)- ( 13 ) 



We refer to |L9| for the expression of the exponent s. In the numerical analysis below 
exp(— s) is set to unity as pj < 1/6. 

The single logarithms in the various subprocesses are summed using the RG equations 



H r (M w ,p) 
H(z,t,fi) 



exp 
exp 
exp 



11 d ^ i i-w 

1/6 /! 3 

f dp, . . 

— jsM/j,)) 

1/6 fl 

^ dp, 



exp 



/M), 

H r (Mw) , 

/ — E7j,(««(^)) +7s(a s (//)) +7//,.(«s(/ i ))] 



H(z,t),(U) 



JM W /i 

rt dp ! 
//- £ , 

7j. = —2a s /Tr, 75 = —CF<y s /n and 7^ r being the anomalous dimensions of the jet function, 
the soft function and the harder function, respectively. Jj(b), f(z,b) and H r (Mw) on 
the right-hand sides are the initial conditions of the RG evolution, and do not contain 
large logarithms. Therefore, Jj(b) are approximated by the tree-level delta functions just 
like Ji(b), and H r (Mw) takes its lowest-order expression H r (Mw) = 1- The B meson 
distribution function f(z, b) is the same as that employed in the semileptonic decays. The 
scale t should be chosen as the typical scale of the hard subprocess in order to minimize the 
involved large logarithms. A natural choice is the multiple of the maximal relevant scales, 

1 



t — K max 



Pj 



■ b 



(15) 



with k an adjustable parameter. Though other scales besides t still lurk in the hard part, 
the logarithms they generate are relatively small. Since t depends on the quark kinematics, 



it reflects the details of the meson dynamics, t remains a hard scale as long as the large 
b region is Sudakov suppressed by exp(— s), and the perturbative calculation of the initial 
condition H(z,t) is reliable. Hence, H assumes its lowest-order expression in Eq. (|]). 

Combining the exponents in Eq. (|14]) , the QCD evolution from Mw to t, governed by 
7# r , is identified as the Wilson coefficients. Conventionally, the summation of the large 
logarithms ln(Mjy//i) is performed in terms of the effective Hamiltonian with W boson 
integrated out. The short- distance physics related to the W boson are incorporated into 
Wilson coefficients. For nonleptonic decays, the relevant effective Hamiltonian (ignoring 
the penguin operators for the moment) is 

4G F _ 



H. 



eff 



V2 



V ch V: s [cM0 1 + cM0 2 } 



(16) 



with the four-fermion current- current operators 



Ox = (s L 'y ll b L )(c L 'fc L ) , 2 = (c L 7 M 6 L )(s L 7 M c L ) 



(17) 



0\ 2 require renormalization and thus the Wilson coefficients ci^in) depend on the renor- 
malization scale //. It is simpler to work with the operators 0± = \{0 2 ± 0\) and their 
corresponding coefficients c± = c 2 (/i) ± Ci(/x), because they are multiplicatively renormal- 
ized. 

At last, the nonleptonic decay rates are written as 



r 



Ml 



dx dy dyo dz db b 
Jo 



N c + 1 



c 2 + (t) + 



N r -1 



c 2 -{t) 



■.Hiz^HJiityHjjityf&b) exp 



^s(pj,6) - s JS (t,b) 



with 



Jl/b fX 



(19) 



The expressions of the Wilson coefficients c± up to next-to-leading order are referred to 
|[20|| . Hence, the three-scale factorization theorem contains the two-stage evolutions: from 
My/ to t due to the summation of the logarithms \n(My//t) organized into the Wilson 
Coefficients, and from t to 1/6 due to the summation of the logarithms ln(£&) into the 
exponential exp(— sjs)- 



8 



We compute the B meson decay rates of the various modes for different choices of the 
hard scale parameter k. The kinematic inputs are Mb = 5.28 GeV, = 4.6 GeV, which 
is determined by the first moment A/ Mb = 1 — nib/ Mb of the distribution function in 
Eq. (|i~0D , Mb = 1.87 GeV, and m c = 1.6 GeV. Our predictions for I?sl and n c are not 
sensitive to the variation of m c . The difference between the values for m c = 1.6 GeV and for 
m c = 1.5 GeV is less than 5%. For the next-to-leading-order Wilson coefficients, we adopt 
the HV scheme. The difference between the HV and NDR schemes has also been examined, 
and found to be less than 5%. In Table 1.1 we list the results in the obvious notations r, 
with the partial decay rates normalized to 5(5 ^ X c cv). This table is compared to the 
HQEFT predictions in Table 1.2 [|J. On the experimental side, the CLEO group reports 



B SL = 10.19 ± 0.37%, n c = 1.12 ± 0.05% and recently B{b -> ccs) = 21.9 ± 3.7% M 



which are only marginally consistent with the LEP measurements B^l = 11.12 ± 0.20% 
and n c = 1.20 ± 0.07%. 

Our predictions show a tendency that the experimental data can be accommodated if 
the hard scale is chosen properly. Unlike the next-leading-order HQEFT calculation, we can 
lower the semileptonic branching ratio into the range of experiment data with the charm 
counting almost remaining fixed. The value of B(b ccs) in Table 1.1 is consistent with 
the CLEO measurement. The b — > cud mode gains larger enhancement than the b — > ccs 
mode in our approach, because of the larger available phase space for the QCD evolutions. 
Therefore, we predict a larger r u d than that from HQEFT. On the contrary, the calculation 
of Bagan et al. shows that the increase of the mode b — > ccs from 0(a s ) QCD corrections 
dominates 0. Though there is an experimantal argument for r ud to be about 4.0 |22 |, the 
scene is still problematic pifl . Indeed, if r ud is about 4.0 and B{b — > ccs) is 21.9%, which is 
determined more accurately than r u d by CLEO, the branching ratios of the various modes 
will not add up to unity To render the total branching ratios equal to 100%, r u d would 
be 5.2 ± 0.6, consistent with our predictions. Our results of the B meson lifetime T, whose 
errors come from the uncertainty of \V c b\ = 0.041 ±0.005, are also consistent with the data 
(1.62 ±0.06) x 10~ 12 sec. 

Note that our calculation leads to r TV = 0.378, which is much larger than the HQEFT 
result 0.23. A simple investigation indicates that r TV depends on the B meson distrubution 
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function. By varying the distribution function, a smaller r TV can be obtained. However, we 
stress that the distribution function employed here is determined from the decay spectrum 
of B — > Xg'-y. Therefore, we urge experimental measurements to decide whether the PQCD 
approach or the conventional HQEFT approach gives better results. 
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Table 1.1 PQCD predictions of the various decay modes with the notations r TV = B(B 
X c tV)/B{B -> X c eu), r ud = B(B -> X C )/B(B -> X c ez>), and r cs = £(£ -> X c5 )/B(B 
X c eV). 



k r TU r ud r cs B c£s B SL n c T (1(T 12 sec) 

1.0 0.378 4.05 1.55 19.41 12.52 1.195 1.50 ±0.37 

1.5 0.378 4.76 1.74 19.59 11.26 1.196 1.35 ±0.33 

2.0 0.378 5.44 1.88 19.36 10.30 1.194 1.24 ±0.30 



Table 1.2 HQEFT predictions of the various decay modes. 



n/m b 




r u d 








n c 


1.0 


0.22 


4.21 


1.89 


22.68 


12.0 


1.20 


0.5 


0.23 


4.75 


2.20 


23.98 


10.9 


1.21 
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